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ABSTRACT 

The problem of estimating the distribution of detection ranges, e.g., of a 
sonar or radir. in which there is no information on the total number of detection 
opportunities or the target tracks, has been treated in OEG lRM-7, (reference 
(1)). A more detailed discussion of the assumptions and derivations is presented 
here; the result differs slightly from that given in lRM-7. 
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We are given an observer with a detection device, and a collection of 
numbers r. , :-i.N, representing the ranges at which targets of a 

certain type have been detected. It is assumed that the detection device has 
a range r, called tlie detection range, which varies randomly as time passes, 
but that it is essentially constant during the passage of any one target. A 
target is assumed to be detected, at range r, if and only if it passes within a 
range r of the detection device. The path of each target is a straight line. 

Only one target at a time passes by. 

It is known that there are cases in which targets pass directly by the 
detector and go undetected. We may think of the detector as "not working" in 
such .cases. It is easy to include this situation in our model; we simply define 
the detection range r to be r=-l when the device Is not working. U i is the 

probability that the detection device is not working, then, of course, r = -1 

with probability i' . We cannot estimate Y in the situation described here, 
o o 

The main thing is the assumption that the detection range against the given 
type of target is a random variable associated with the detection device; it does 
not depend on the individual target, or indeed, on the presence or absence of a 
target. During the total time of the "experiment, " N targets were detected, at 
ranges r, , .... ; we don’t know how many targets passed by that were not 

detected. 

We shall denote by o the lateral range of any target (detected or not); this 
is the shortest distance from the observer to the straight line path traversed by 
the target, l.e,, the "distance of closest approach." From the definition of the 
detection range, we see that a target Is detected if and only if os ps r. Both 
0 and r are random variables; o is associated with the target, and r with the 
detection device. It follows that p areT are independent random variables. 

In the following, r^ denotes an arbitrary, but fixed, positive number. 

We are inreresied in two functions: 














I 


f(r^)^Pr(r^r^). 


and 


( 1 ) 

l 


L(r ) = Pr(a target is detected, given that its lateral range 

° is r ). 

□ 


( 2 ) 
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The firsi function depends only on the doteecion device, but we can- 

estimate it only through observed detections, i.e., targets arc needed. 

The probability space that is involved here is evidently the (;. , r) plane, 
p ^ o , r 2; -1. The observed ra.nges r^, .. of the detected targets are 

independent samples of a certain probability distribution. This distribution is 
nothing more nor less chan Che distribution of detection ranges (of the detection 
device) given a (target) detection, i.e., the distribution of r given os 2 5. r, 
if we recall that a target is detected if and only if o s: p s r. 

• . Regarding (2), we can write 


L(r^) = Pr(o s c ^ > 0 ' = 


= Pr(o i r^ i r I 0 = r^) . 

But in the last e-xpression, the two events r s r and q = r are independent. 
Hence . • “ * ° 


L(r^) = Pr(r^ s r) = f(r^) . 


» • 


by (1). This appears 


afs in 


as Theorem 1. 


So all we need nqw is a reasonable estimate for Since we have 

plenty of .samples-of the distribution of r given detection, we can write 


Pr(r sr<r +.Ar I detection)^ , 
0 0 o N . 


in which An is the number of observed detections with ranges in the interval 

(r ,r-t-Ar),r > 0 , and N is the total number of detections. But 
000.0 


Pr(r sr<r +Ar Idetection) 

*0 o o 

■=Pr(r si'<r + Ar l o s n ^.r) 

' .0 o o 

' Pr(r :£ r < r ’+ Ar and o ^ r) . 
o o o 
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Let a = Pj^(o * r) = Pr(decectiQn). 

Now the numerator of (5) differs only negligibly from 

Pr(r^ s r < r^ + Ar^ and o ) 

» 

if Ar^ is small, and provided that q has a continuous distribution. 

Since -r sr<r +Ar | and Xq £ r !• are independent events, (4) and (5) 
o o o , L o J 

now yield 

Pr(r £ r < r + Ar). PcCq £ r ) An ^ 

O O O 1^2 V t • 



or 


Pr(r £r<r +Ar)as 

' O 0 0 


N 


An 

Mpi r ) 


( 6 ) 


. Now let us suppose that q is uniformly distributed in an interval 
o s 0 ^ Then (o) becomes 

Pr(r £r<r +Ar)=-^ 

' o o o N 


if r < R . 
o 


(7) 


-V 7 An if r a R . 
N o 


la aa uapublished working paper, J. Neueudorffer, who first considered 
this problem, groups his observed ranges into range bands. While his I 

estimate is basically good. It contains an artificial bias that is introduced by 
the choice of end points for his range bands. 

Here we shall simply take as our estimate for ~ or rather 1 - 

which is nondecreasing that discrete distribution which is determined by 

the right side of (7). That equation will be satisfied if we take the jumps in I 

1 - f(r) — we drop the subscript from r — to be located hi r = -I and 
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r = , __ , witli niagaituJe.s ^ . ... ir 


ia which 


= Pr(r = -1), unk nown, 


Y - -3- 
i ■ N 


if r. > R . 


By summing for i = 1.N, we obtain 


1 - y = - ° [ y -5. 


rj s R i 


in which M is the number of observed r. larger than R. From (9). we can 
estimate o/N; ‘ 


a/N =- 


z: 

1 - s R •» 


» • 


and (8) can be written 


ifo =Wr = -l), 


S'. ■ -^--T 


‘ » R 4 . 

— + M 


®i “ V o 

iL. + 

Tj S R ’"j 


if r . s R 
1 


if r. > R , 


B 
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and finally, _ ^_ 

f(r)= E 5^ = L(r). ■ j 

r. i r ' ' 

1 

This estimate, for rio, contains the unknown Y in the constant 
factor 1 - J , which does not affect the shape of the curve. 
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